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We study the non-equilibrium condensation process in a holographic superconductor with nonlin¬ 
ear corrections to the U(l) gauge field. We start with an asymptotic Anti-de-Sitter(AdS) black hole 
against a complex scalar perturbation at the initial time, and solve the dynamics of the gravitational 
systems in the bulk. When the black hole temperature T is smaller than a critical value T c , the 
scalar perturbation grows exponentially till saturation, the final state of spacetime approaches to 
a hairy black hole. In the bulk theory, we find the clue of the influence of nonlinear corrections in 
the gauge filed on the process of the scalar field condensation. We show that the bulk dynamics 
in the non-equilibrium process is completely consistent with the observations on the boundary or¬ 
der parameter. Furthermore we examine the time evolution of horizons in the bulk non-equilibrium 
transformation process from the bald AdS black hole to the hairy AdS black hole. Both the evolution 
of apparent and event horizons show that the original AdS black hole configuration requires more 
time to finish the transformation to become a hairy black hole if there is nonlinear correction to the 
electromagnetic field. We generalize our non-equilibrium discussions to the holographic entangle¬ 
ment entropy and find that the holographic entanglement entropy can give us further understanding 
of the influence of the nonlinearity in the gauge field on the scalar condensation. 

PACS numbers: 


I. INTRODUCTION 


The AdS/CFT correspondence has been proved as a unique approach to study the strongly coupled 


field theories 1 3]. 


Recently, there appeared a new direction in applying this correspondence to shed 
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some lights in condensed matter physics, with the hope that the strongly coupled systems can be 
computationally tractable and conceptually more transparent (for a review, see (4-6] and references 
therein). The simplest example of the duality between the superconductor and gravity theory was first 
formulated in [7j. In the gravity side one has Einstein-Maxwell-charged scalar theory with negative 
cosmological constant, which usually has static Reissner-Nordstrom-AdS (RN-AdS) black hole solutions 
with vanishing scalar hair. However in the low temperature, because of the coupling between the 
Maxwell held and the charged scalar Held, the RN-AdS black hole is expected to become unstable and 
finally the black hole will have the scalar hair breaking the U(l)-gauge symmetry spontaneously. The 
RN-AdS and hairy black holes are remarkably identified as the normal and superconducting phases 
of a superconductor realized in the dual Held theory on the AdS boundary. This striking connection 
between condensed matter and gravitational physics is expected to be a hopeful approach to disclose the 
properties of strongly correlated electron systems. Most available studies on this topic were concentrated 
on the static AdS black hole backgrounds dual to the superconductors reflected by the dual operators 


on the AdS boundary with equilibrium or nearly equilibrium dynamics 8]- 


33]. 


As gravitational theorists, one of our motivations of seeking a dual description of superconductivity, 
which is below our usual energy bandwidth, is to understand more on the black hole physics, for example 
how the scalar hair attaches to the bald AdS black hole to make the black hole become hairy step by 
step, what conditions stimulate or prevent the condensation of the scalar hair onto the AdS black hole, 
what characteristics of the transformation between the AdS bald black hole and the hairy one are, 
what spacetime properties are for the developed hairy AdS black holes, etc. To answer these questions, 
we cannot just look at the dual operator on the AdS boundary in an equilibrium system, we need 
to concentrate on the non-equilibrium dynamics in the bulk gravitational configurations to disclose 
the process of addressing the scalar hair onto the bald AdS black hole. In 34;], the non-equilibrium 
process in the bulk Einstein-Maxwell-scalar background has been studied. It was observed that during 
the time evolution, the initial scalar perturbation in the bulk spacetime first moves towards the AdS 
boundary, it is then bounced back by the AdS wall, grows and moves towards the bald black hole. 
Subsequently this perturbation falls onto the bald black hole and accumulates there and finally makes 
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the bald black hole become hairy. The time evolutions of the event and apparent horizons of the bulk 
black hole in the AdS bulk were disclosed in 34J. Employing the AdS /CFT duality, the authors related 
the observed properties in the AdS bulk to the boundary CFT, and further they displayed the time 
dependence of the superconducting order parameter in the holographic superconductor. The attempt 
in the non-equilibrium approach has been generalized to the study of the dynamical phase diagram 
corresponding to quantum quenches of source field where the bulk black hole perturbation spectrum 
has been investigated in relating to the distinct behaviors on the boundary CFT 35]. More extensions 


of the non-equilibrium physics in describing the transformation of an unstable anisotropic black hole 
without scalar hair into a stable hairy black hole, which relates to the holographic superconductor with 
anisotropy, have been done in 36]. Further works on the far-from-equilibrium dynamics of holographic 


superconductor including periodic driven system [37], quantum quench of a superconducting AdS soliton 

39, 4t] 


38j , holographic thermalization 


and the formation of topological defects 


41 


42] have also been 


extensively studied. 

Motivated by the application of the Mermin-Wagner theorem to the holographic superconductors, 
there have been a lot of interest in generalizing the Einstein gravity background in the holographic 
superconductor to the gravity including the curvature correction. By examining the char ged scalar field 


together with a Maxwell field in the Gauss-Bonnet-AdS black hole background 21,|22j, 43H52|. it has 
been observed that the higher curvature correction makes the holographic superconductor harder to 
form. Besides the curvature correction to the gravity, it is also interesting to investigate the high-order 
correction related to the gauge field. In the study of the ratio of shear viscosity to entropy density, 
it was found that the higher derivative correction to the gravity has different effect from that of the 


higher derivative correction to the gauge field 


53]. This motivates people to study the effect of higher 


order corrected Maxwell field on the scalar condensation and compare with the effect brought by the 
correction in the gravity 


54|. 


Since Heisenberg and Euler 55| noted that quantum electrodynamics predicts that the electromagnetic 
field behaves nonlinearly through the presence of virtual charged particles, the nonlinear electrodynamics 
has been studied a lot for many years [55-61]. In the low-energy limit of heterotic string theory, the 
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higher-order correction term appears also in the Maxwell gauge field 


56J. Born-Infeld electrodynamics 


57| is the only possible nonlinear version of electrodynamics which is invariant under electromagnetic 


duality transformations 


58j . The nonlinear electromagnetic fields in logarithmic form appears in the 


description of vacuum polarization effects. The form was obtained as exact 1-loop corrections for 


electrons in an uniform electromagnetic field background by Euler and Heisenberg 


55j. In 


621, the 


author discussed a simple example of a Born-Infeld-like Lagrangian with a logarithmic term which can 
be added as a correction to the original Born-Infeld one. Replacing the Maxwell field with the nonlinear 
electromagnetic field, it was disclosed from the behavior of the boundary operator in the CFT that 
the nonlinearity in Born-Infeld and the logarithmic form fields can make holographic superconductor 
difficult to form [19L 63]. In the corresponding bulk spacetime, it was observed that the nonlinearity of 


the electromagnetic field can hinder the growing mode to appear in the scalar perturbation 64], which is 
consistent with the phenomenon on the boundary CFT. Other works on holographic dual models in the 
presence of nonlinear corrections to the Maxwell electrodynamics can also be found in [65-67]. These 


discussions were concentrated on the equilibrium state both in the bulk and boundary. To understand 
more on how these two kinds of nonlinear corrections in the electromagnetic fields influence the scalar 
condensation in the holographic superconductor, we can not only concentrate on the equilibrium state 
but also study the non-equilibrium condensation process in holographic superconductor with nonlinear 
electrodynamics. This can help us see closely how the nonlinearity in the electrodynamics influences 
the condensation of the scalar hair in the process of the black hole evolution and transformation, it also 
can aid us to see how the nonlinearity affect the phase transition in the boundary system . 

In this paper we will concentrate more on the non-equilibrium bulk spacetime. The organization of 
the paper is the following: In Section II, we will introduce the holographic superconductor model with 
nonlinear electrodynamics. We will give equations of motion and the initial and boundary conditions to 
solve them. In Section III, we will show the numerical results in the non-equilibrium bulk system. We 
will present that the condensation process of the scalar field depends on the strength of the nonlinearity 
in the electromagnetic field. We will describe the time evolution of the event and apparent horizons of the 
black hole in the bulk. Furthermore, we will display the time evolution of the holographic entanglement 
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entropy and its dependence on the nonlinearity of the electromagnetic field. We will summarize and 
discuss our results in the last section. 


II. MODEL 


The action describing a gauge field coupled with a scalar field with a negative cosmological constant 
in the 4D spacetime reads 


S= d A Xy/—g 


7? + -| + L(F ) - IVS' - iqA^ I 2 - m 2 |4'| 2 
V 


(1) 


where l is the AdS radius which is set to unity hereafter, q is the coupling factor between the gauge field 
and the complex scalar field, m is the mass of the scalar field. L{F ) is the Lagrangian density of the 
nonlinear electromagnetic field. In this work, we will investigate the nonlinear electromagnetic fields in 
two different forms, the Born-Infeld form (BINE) and the Logarithmic form (LNE), 

1 


L(F) = 


b 2 ^-V 1 + ^ 2F2 J’ bine, 

- Ihnfl+ ±b 2 F 2 Y LNE, 


where F 2 = F 1 ' 1 'F IIU . and b is the nonlinear parameter. As b —> 0, both kinds of the nonlinear electro¬ 
magnetic fields reduce to the Maxwell field, L(F) —> —F 2 /4, and the model action dTJ) reduces to that 


in [341. 

We take the ansatz of the background as the ingoing Eddingthon-Finkelstein coordinates, 
ds 2 = — — [f(v, z)dv 2 + 2 dvdz\ + (j) 2 (v, z){dx 2 + dy 2 ). 


( 2 ) 


where v is the ingoing Eddington-Finkelstein time coordinate, z is the radial direction in the bulk and 
the AdS boundary is located at z = 0. With the axial gauge, the ansatz of the gauge field and scalar 
field can be written as, 


A = (a(v,z), 0,0,0), 


= il)(v, z). 


( 3 ) 


There are still the following residual gauge symmetries for the anstaz which will be fixed by the boundary 
conditions 

>^-+g(v), a ^ a + d v 9(v), if) —> 
z z 


( 4 ) 
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From the action and the ansatz we can derive the equations of motion as follows, 


(<W + ^2 (6 - m 2 \ip\ 2 + L(F) - l -z A F 2 a' 2 ) = 0, 

(5) 

Dip' + ~jpDip + + Tjiqat'ip + ^ip = 0, 

(6) 

(z 2 (4)')' - z 2 a ,2 F 2 + 4^D<P - {iP'*DiP + iP'DiP*) = 0, 

(7) 

2 z 2 F 2 {Da)' + z i (z- 2 fF 2 )'a' + 4 z 2 F 2 ^a' 

+2z 2 DF 2 a - z 2 fF' 2 a - 2iq{ipDip* - ip*Dip) = 0, 

(8) 

D^+l^'DcP + ^DiPl^O, 

(9) 

( Z 2 <P'Y + y^lV’T = 0, 

(10) 

(z 2 F 2 a'Y + 2z 2 F 2 a'— + iq(ip'*ip - ip*ip') = 0 

<P 

(11) 


where the prime denotes the derivative with respect to z, i.e. 1 = d z , and the derivative operator D is 
defined as 


Dtp = d v tp-fd z tp/ 2, 

Dip = d v ip — fd z ip/2 — iqaip, 
Da = d v a — fd z aj 2, 


D 2 <t> = d v (Dc/>)-fd z (DcP)/ 2, 
Dip* = d v ip* — fd z ip*/2 + iqaip*, 
DF 2 = 8 V F 2 - fd z F 2 / 2. 


The function F 2 is defined as F 2 = (1 + b 2 F 2 / 2) -1 / 2 for BINE model, while for the LNE model it reads 
F 2 = (1 + b 2 F 2 / 8) _1 . In the limit b —> 0, F 2 —> 1, the system and the equations of motion go back to 
those described in [34]. In our numerical computation, we will mainly consider the first four equations, 


Eqs. (0 —( 0 ), describing the evolutions of the variables /, (p , ot, and ip. We will treat the remaining three 
equations as constraints. 


A. the normal phase 

While the scalar field is null, the solution to the equations of motion is a static AdS black hole solution. 
As the nonlinear parameter b approaches zero, £> —>- 0, the nonlinear electromagnetic field goes back to 
the Maxwell field, the background is nothing but the Reisser-Nordstrdm-AclS black hole, the solution is 
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given by 


= —, ■0 = 0, a = Qz, f = 1 — 2Mz 3 + -Q 2 z 4 . 

z 4 


( 12 ) 


While b ^ 0, one can solve the equations of motion directly, the solutions are given as follows: 
( i ) In the case of BINE 


<t> = ~, V> = 0, a = Qz 2 F 1 [j, i ^,-b 2 Q 2 z 4 ] 

Z 4 Z 4 

/ = 1 - 2 Mz 3 + ^2 (1 - \A + ^ 2 Q 2 2 4 ) + 2 -Fi[~, -& 2 QV]. 


(13) 


[ii) In the case of LNE 


= —, 0 = 0, a = 

z 


3i4? ( i~ ' /1 + !,2Q2z4) + 5 QZ2Fl[ i’5'r~ i,2QV1 ' 


/ = W2M?+- (1-01 + 6 2 Q 2 0) 


+ 9 2Fl[ 4’2’4’ b2Q2z ^ 3 0 Ln( ^ 6 2 Q 2 0 


1 0(01 + & 2 Q 2 0 — 1), 

3b 2 


(14) 


where 2 -P 1 represents a hypergeometric function. Q and M are proportional to the charge and the mass 
in the field theory, respectively. The spacetime described in Eq. (fl3l) is nothing but the Born-Infeld-AdS 
black hole 6f| and the spacetime presented in Eq. (fill) is the same as that in 
For these static black holes, the Hawking temperature is given by 

1 df. 


621. 


T= — 


4-7T dz 


(15) 


Zh is the radius of black hole horizon defined by f(zh) = 0. These asymptotic AdS black holes are 
unstable when the temperature of the black hole T is smaller than a certain critical value T c . The 17(1) 
gauge symmetry will be spontaneously broken due to the condensation of the complex scalar field. 












B. Boundary and initial conditions 


We set the mass of the complex field as m 2 = — 21 2 . The expansion of the variables near AdS 
boundary behave as 

f(v,z) = 1 +fi{v)z +f 2 {v)z 2 + f 3 {v)z 3 + ..., 

fav,z) = 1/z + fa(v) + <t>l(v)z + fa{v)z 2 + fa(v)z 3 + ... , 

fav,z) = fa(v)z + fa{v)z 2 + +fafa)z 3 . . . , 

a(v,z) = ao{v) + a\(v)z + a 2 (v)z 2 + a 3 (v)z 3 - (16) 

With the residual symmetries Eq. 0, we can choose both fa{v) and ao(v) to be zero. Substituting Eq. 
m into Eqs. (151151) . solving them we can obtain the coefficients as the functions of fa{y) and fa(v), 

fi =0, f 2 = — fa = -^#i| 2 ; 

fa = ~(fafa + fafa), fa = ~^\fa \ 2 ~ ^l^il 4 - \{fafa + fafa), 

a .2 = ~iq{~fafa+fafa), fa = + fa + rV’iV’i! ( 17 ) 

together with other two formulas, 

h = + fa$* ~ fafa ~fafa), d 3 = -iq(fafa - fafa + fafa - fafa ). (18) 

where the overdot denotes the derivative with respect to v ,’ = d/dv. a 3 and f 3 are related to the charge 
and mass density of the system, respectively. To solve the system numerically, we employ the Chebyshev 
Pseoduspectral representation in ^-direction, and an implicit finite difference scheme in redirection. We 
found that 20 points in ^-direction are sufficient to assure reliable results. The data presented in this 
paper is for 30 grids. Convergence testing shows the errors will be less than one percent in the plotted 
quantities. 

In our numerical process, we consider fafa) as the source of the boundary field and set it to be null 
at the AdS boundary z = 0. In Eq. m, it is interesting to see that if one chooses fa = 0, f 3 and a 3 
are conserved quantities in the system, which means that the black hole charge and mass are constants 
in the phase transition process. On the initial moment we specify the black hole M and Q as fixed 
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parameters since they are conserved quantities in our setting. Without loss of generality, using the 
following scaling symmetry, 


(v, z, x, y) —> ( kv , kz, kx, ky), 

F—y F, 4>—y <t>/k, a—> a/k, ip — yip, 

M — yM/k 3 , Q — yQ/k 2 , T—y T/k. 


we can set the initial black hole horizons radius to unity. With /(1) = 0 we can get 
(i) In the case of BINE 

M = ) [1 + eP (1 - j, -Wit 


(ii ) In the case of LNE 


M = 2 [ 1 + 9^ (1 ' x/1 + W) 


4Q 2 z 


4 F' 1 
2-^1 


1 5 , W1 _ 1 Lnf 2( ' /TT ^- 1) )l 

4’2’4’ ^ k 2 Q 2 ■' 


thus only the black hole charge Q is left as the open initial variable to be assigned. As 
condition, the background is given by the normal charged black holes, described by Eqs. 
null scalar hair, along with a scalar perturbation in the following form, 

ip(z) = z 2 Aex p 


(z - Zm) 
2S 2 


for the initial 
(11211141) . with 


(19) 


with A = 0.01 ,5 = 0.05 and z m = 0.3. 

III. NUMERICAL RESULTS 


In this section, we will show the numerical results for the non-equilibrium process for the holographic 
superconductor with nonlinear electrodynamics. We will first exhibit the dynamics of the scalar field, 
including its condensation in the bulk and the evolution of the order parameter on the AdS boundary. 
Besides, we will disclose the evolution of both the event and apparent horizons. Furthermore, we will give 
the evolution of the holographic entanglement entropy in this time-dependent background. Hereafter, 
we take q = 4 in our numerical calculation. 
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A. The Condensation of Scalar Field 


To calculate the critical temperature T c , we use the shooting method the same as those in [161.121.151]. 


With this method, we can reproduce the phase transition temperature in 


34j. The critical temperatures 


of black holes with nonlinear electromagnetic field in the forms BINE and LNE are listed in Table. I. We 
see that for both nonlinear electromagnetic field models, the critical temperature T c decreases with the 
increase of the nonlinearity. This tells us that the nonlinear correction to the Maxwell gauge field hinders 
the formation of the hairy black hole, which also means nonlinearity make the boundary system harder 


to reach the superconducting phase. The qualitative behavior of T c is consistent with that disclosed in 


the static black hole backgrounds in [63 ], 


b 

0 

0.6 

1.2 

1.5 

1.8 

T c 

BINE 

0.2195 

0.2112 

0.1879 

0.1720 

0.1543 

LNE 

0.2195 

0.2154 

0.2044 

0.1971 

0.1889 


TABLE I: The critical temperature T c as a function of the nonlinear parameter b for both black holes with BINE 
and LNE with the backreacting parameter q = 4. 


In FigQ] starting from the initial black hole temperature T ~ 0.5T C , we depict the time evolution 
of the scalar field q\ip{v, z)\ as functions of vT c and zT c for different strengths of the nonlinearity in 
the electromagnetic field. The general feature of the condensation process is similar to that reported in 
34]. The initial Gaussian perturbation runs to the AdS boundary, then is reflected into the bulk and 
keeps growing exponentially until the saturation is reached. In the left panel, the parameter b is fixed 
as b = 0. We can see that the scalar field keeps on increasing exponentially within the time interval 
0 < vT c < 5, but starting from vT c ~ 5 the scalar field becomes saturated at the horizon. The saturation 
of the scalar hair at the black hole horizon marks that the gravitational system enters into the phase 
of a static hairy black hole state. At the saturation moment, the scalar field q\ip{v, z)\ near the black 
hole horizon reaches q\ip{v, z)\ ~ 1.3. In the middle panel, we plot the figure with the parameter b = 1.8 
for the BINE model, while in the right panel we exhibit the result for the LNE model with the same 
parameter b = 1.8. Comparing the left panel with the middle and right panels, we see that with the 
increase of the nonlinearity in the electromagnetic field, the scalar field needs longer period to reach the 
saturation, vT c ~ 6 in the middle panel. Besides we see that for the electromagnetic field with stronger 
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nonlinearity, the final stable hairy black hole has more scalar hair, higher value of q\ip(v, z)\ (for example 
q\ip(v, z)\ ~ 1.8 in the middle panel). 




FIG. 1: The scalar field \ip(v,z)\ as functions of vT c and zT c , it grows exponentially till saturation. In the left 
panel the parameter 6 = 0. The middle panel is for a system with BINE, and the right panel is for a system 
with LNE, b is fixed as b = 1.8. 


B. Evolution of the order parameter 


In this subsection, we describe the dynamics of the order parameter of the boundary theory. Following 


m 3 , the order parameter is defined as 


( 0- 2 {y)) = V2Mv). 


( 20 ) 


In FigO which displays the raw data, we show the numerical results with initial black hole temperature 
as T ~ 0.5 T c . It depicts the time-dependence of the order parameter (q\(02(v))\) 1 ^ 2 /T c with b = 
0,0.6,1.2,1.5 and 1.8 respectively, for both BINE and LNE models. From this figure, we see that no 
mater what the value of b is, given enough time the order parameter will reach some certain values 
and will not increase any more, which implies that the system in the boundary CFT evolves into an 
equilibrium phase, i.e. the superconducting phase. The behaviors displayed in Fig. [2] on the order 
parameters of the CFT on the boundary support the properties of the scalar hair evolutions disclosed 
in the bulk, for the model with stronger nonlinearity in the electromagnetic field, it costs more time 
for the order parameter to approach the saturation, the boundary is more difficult to arrive at the 
superconducting phase state. Fig[3] shows the curves normalized by T c . With respect to Fig|2] along 
the vertical direction the effect of T c is to increase the height of these curves, but it doesn’t differ the 
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order of these curves from bottom to top, also the properties of the curves. In horizontal direction, the 
effect of T c is to decrease the time period of non-equilibrium process in the boundary system, but the 
time interval with a bigger b is still longer than that with a smaller b. 




FIG. 2: This figure depicts the evolution of scalar operators in the AdS boundary as a function of v/L with 
initial temperature T ~ 0.5T C and different values of the nonlinear parameter b. The left panel is for a model 
with BINE, while the right panel is for LNE. 




vT t . 


vT c 


FIG. 3: This figure depicts the evolution of scalar operators in the AdS boundary as a function of vT c with 
initial temperature T ~ 0.5T C and different values of the nonlinear parameter b. The left panel is for a model 
with BINE, while the right panel is for LNE. 


C. Evolution of Horizons 

To further disclose the properties reported in the previous two subsections, we extend our discussion 
here to the evolution of black hole horizons. 

The apparent horizon is defined as the outer boundary of any connected trapped region. This definition 
makes the apparent horizon a marginally outer trapped surface on which the expansion of outgoing null 
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geodesics vanishes. Following this definition, the apparent horizon can be obtained by 


D4>(v,z A h(v)) = 0. 


(21) 


The event horizon is spanned by the outgoing light rays which neither hit the black hole singularity nor 
reach future null infinity. An obvious approach to locate the event horizon is to examine null geodesics 
34j by solving the equation 


Zeh(v) = z E h{v ))• 


( 22 ) 


As we get the zah{v ) and zeh{v), following 
horizons as 


34l | we can calculate the area of apparent and event 


Area(apparent/event horizon) = cj>(v , zah/zh{v)) 2 . 


1.15 


a 1.10 
< 

1.05 

1.00 


FIG. 4: This figure depicts the area of apparent and event horizons as functions of v/L with initial temperature 
T ~ 0.5T C for various nonlinear parameter b. The left plot is for a holographic superconductor with BINE while 
the right one is for a model with LNE. Solid lines correspond to event horizons while dashed lines apparent 
horizons. 

Fig (4| displays the time evolution of the area L 2 Area of the apparent and event horizons as functions 
of v/L in a gravitational system with nonlinear electromagnetic fields described by BINE and LNE 
models, respectively, where we choose the initial black hole temperature T ~ 0.5T C . The initial black 
holes’ horizon radius are set to unity, thus for an initial static AdS black hole the initial area of event 
horizon, which also equals the black hole entropy, coincides with that of apparent horizon. Furthermore, 
we can clearly see the transformation. The original stable bald AdS black hole with constant entropy 
evolves into the final stable hairy black hole with higher constant entropy. This transformation satisfies 
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vTc vTc 

FIG. 5: This figure depicts the area of apparent and event horizons as a function of vT c with initial temperature 
T — 0.5 T c for various nonlinear parameter b for a holographic superconductor with BINE. Solid lines correspond 
to event horizons while dashed lines apparent horizons. 



vTc vTc 

FIG. 6: This figure depicts the area of apparent and event horizons as a function of vT c with initial temperature 
T — 0.5 T c for various nonlinear parameter b for a holographic superconductor with LNE. Solid lines correspond 
to event horizons while dashed lines apparent horizons. 


the second law of thermodynamics. The increasing of the area of the horizon exhibits the non-equilibrium 
process in the transformation. From both the two plots, it is clear that when the nonlinearity of the 
electromagnetic filed becomes stronger, the horizons need longer time to settle down to a constant value, 
which means that it needs more time to develop a stable hairy black hole. Further we see that the final 
stable state develops higher entropy in the gravitational system with stronger nonlinear electromagnetic 
field. 

Figs (5] and [5] displays the normalized data. The left plots present the evolution of the area Area/T c 2 . 
At vT c = 0, for the initial black holes with different 6, the entropy no longer takes the same value, this 
is clearly because the impact of increasing nonlinear parameter b is that the critical temperature of the 
system is lowered. To see more closely on the influence of the nonlinearity of the electromagnetic filed 
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on the horizons’ evolution, we plot in the right panels the net change of the horizons’ area with respect 
to the value of the initial bald black hole, AArea/T c 2 , in the transformation process. 


D. Holographic Entanglement Entropy 


Finally we extend our discussion to the holographic entanglement entropy. The entanglement entropy 
is considered as an effective measure of how a given quantum system is strongly correlated. It has 
been proven as a powerful tool to describe different phases and the corresponding phase transitions in 
holographic superconductors fl -Q. In this section, we will calculate the holographic entangle entropy 


in the non-equilibrium process and further disclose that the entanglement entropy can be regarded as a 
useful tool in describing the holographic superconductor. 

The entanglement entropy is directly related to the degrees of freedom of a system to keep track 
of them. From [69|, |82|, |83[, we can see in a time-dependent background the entanglement entropy 
associated with a specified region on the boundary in the context of the AdS/CFT correspondence is 
given by the area of a co-dimension two bulk surface with vanishing expansions of null geodesics. We 
will follow this treatment to calculate the entanglement entropy in the time-dependent background. We 
use a strip geometry by considering the length of the strip, L , approaches infinity along the y-direction 
and stretches from x = — h/2 to x = h/2 in the ^-direction, and the width of the strip is h. We express 
v and z by using parametrization v = v(x) and 2 = x(x). In the dual gravity side, this leads to the 
following boundary condition: 


’(h/2) = v(—h/2) = vo, z(h/2) = z(—h/2) = zo, v'(0) = z'(0) = 0, 


(23) 


where zq —■► 0 is the UV cut-off and uo is the boundary time. 

For the reflection symmetry of the background, v(x) = v(—x) and z( x) = z(—x), the resulting 
entanglement entropy S is given by: 

r h / 2 rl T , _ 

4 S/L= / — \Jz 2 <t>(v, z) 2 — f(v, z)v' 2 — 2 v'z 1 , (24) 

J-h/2 z 

where the prime ' denotes the derivative with respect to x, i.e. ' = d x . This integral is independent of 
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x, we have a conserved quantity 

~ ^2^ = z2< t>( V ’ Z ) 2 ~ f( V > Z ) V ' 2 - 2v ' Z: 'i ( 25 ) 

where </>* is short for <j>(v*, z*), v » = u( 0 ) and z* = z( 0 ) are constants. 

From Eq. (1251) . we can derive two equations of motion for v and z, 

fv' 2 — 1 + 2 v'z' + z[cj)(fv' 2 — 2) + v" — v ,2 d z f / 2] + z 3 (f>d z (j)( 2 + z<j>) 

—z 2 [cj) 2 + (j){v' 2 d z f / 2 — v") + (2 v'z' — 1 )d z <f> + 2 v ,2 d v </>\ = 0 (26) 

(1 + z4>){z" + v'z'd z f + v' 2 d v f / 2 ) - 2z 2 d z (z())) 

+zd v (f>[( 1 + 2^>) 2 — 2 'i/ 2 /] + /[u 7/ + cf>(zv" — 2v'z') — 2 zv'{z'd z (j) + v'd v <j>)\ = 0 (27) 

Only two of Eqs. dL m and are independent, we can use two of them to solve v(x) and z(x). 
Finally, we can get the entanglement entropy density S = S/L: 

f h/2 6(v z) 2 

2S(h,vo) = / —7 - dx — Ln(2/zo). (28) 

Jo V* 

In our calculation, we choose the strip width to be 1. To show the dependence of the entanglement 
entropy on the parameters more clearly, we define AS = S(v o) — Sq, where So is the entanglement 
entropy of the strip region in the normal bald AdS state. 



FIG. 7: The holographic entanglement entropy is depicted as a function of Vo/L with the initial tempearture 
T ~ 0.5T C for various values of the nonlinear parameter b in BINE. 

Figs. [7] and [5] depict the raw data of holographic entanglement entropy as a function of the boundary 
time vo/L for various nonlinear factor 6 = 0, 0.6,1.2,1.5 and 1.8 with initial temperature T — 0.5 T c for 
the model of BINE and LNE, respectively. In the left plots, the curves with different 6 intersect with each 
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vo/L vq/L 


FIG. 8: The holographic entanglement entropy is depicted as a function of Vo/L with the initial tempearture 
T ~ 0.5T C for various values of the nonlinear parameter b in LNE. 




v 0 Tc v 0 Tc 


FIG. 9: The holographic entanglement entropy is depicted as a function of voT c with the initial tempearture 
T ~ 0.5T C for various values of the nonlinear parameter b in BINE. 


other, to see clearly the quantity change of entanglement entropy S/L, we plot the net change A S/L in 
the right panels. We see that below the critical temperature when the scalar field starts to condensate, 
the entanglement entropy becomes smaller in the non-equilibrium process. This is consistent with the 
expectation that in the condensation of the scalar field, certain kind of “paring mechanics” triggers the 


decrease of degrees of freedom due to the formation of “Cooper pairs” 78]. When the nonlinearity of 
the electromagnetic field becomes stronger, the non-equilibrium process is longer and the entanglement 
entropy takes longer time to decrease to a constant value. This shows that, fixing T T ^ , in a gravitational 
system with higher nonlinear electromagnetic field, it needs longer period for the system to form the 
“Cooper pairs” to complete the condensate. Fixing T/T c (b) — 0.5, the final value of the entanglement 
entropy for stronger nonlinearity in the electromagnetic field is lower, which shows that for the stronger 
nonlinear electromagnetic field there is a larger reduction in the number of degrees of freedom. 


Figs. [3] and [TO] depict the normalized holographic entanglement entropy as a function of the boundary 
time vqT c . In the left plots, the curves with various b don’t intersect with each other. This is mainly 
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FIG. 10: The holographic entanglement entropy is depicted as a function of voT c with the initial tempearture 
T ~ 0.5T C for various values of the nonlinear parameter b in LNE. 

because increasing parameter b the critical temperature of the system is lowered, as those shown in 
Table. I, thus the impact is that the curves with larger b move upward much more than those with 
smaller b. In the right plots we plot the net change A S/T c . From figs. l9l and fltil we also reach the 
conclusion that, fixing y-yjy — 0.5, in a gravitational system with higher nonlinear electromagnetic field, 
it needs longer period for the system to complete the condensate, and for superconducting phase more 
degrees of freedom have condensed. 



IV. CONCLUSIONS 


In this paper, we have investigated the non-equilibrium condensation process in the holographic super¬ 
conductor with nonlinear electromagnetic fields. We have solved the time evolutions of the gravitational 

2^-3 


system in the bulk configurations. Following 


the asymptotic AdS black hole, which was a static bald AdS black hole at the beginning. When the 
temperature is lowered below the critical value T c , the initial scalar perturbations grow exponentially 
and the original bald AdS black hole spacetime gives way to a hairy black hole step by step. Finally 
the condensation of the scalar hair forces the formation of a new hairy black hole. 

In this non-equilibrium process, we have observed the formation of the scalar hair step by step. 
Especially in the evolutions of the bulk dynamics we have found the clues why the nonlinearity in 
the electromagnetic fields hinders the formation of the hairy black hole. The final configuration of 
the gravitational system indicates that to be stable the black hole needs developing more scalar hair 


in the model with higher nonlinearity electromagnetic field. We also observed that, fixing T 1 { ^ , the 


























19 


nonlinearity makes the condensation process of the scalar field onto black hole longer. In the gravity 
side we also have examined the evolution of horizons in the non-equilibrium transformation process from 
the bald AdS black hole to the AdS hairy one. Fixing T T ^ , both the evolutions of the apparent and 
event horizons showed that with the nonlinear corrections in the electromagnetic field, the original AdS 
black hole configurations require more time to finish the transformation to become hairy black hole. 

As for the superconducting order parameter in the boundary theory, we clarified for the model with 
stronger nonlinearity in the electromagnetic field, it costs more time for the order parameter to approach 
a bigger saturation value, which means the boundary system is harder to arrive at the superconducting 
phase. As we know, for the holographic superconductor the superconducting phase of the boundary 
system corresponds to a stable hairy black hole in bulk. The results obtained in both the bulk and 
boundary sides are consistent. 

Considering the intriguing physical meaning of the holographic entanglement entropy in studying the 
quantum system and its effectiveness in describing the degrees of freedom in a system, we have general¬ 
ized our non-equilibrium study of the problem by examining the holographic entanglement entropy. For 
the stronger nonlinearity of the electromagnetic field, fixing T 1 * 3 [ ^ , the non-equilibrium process is longer 
and the entanglement entropy takes longer time to decrease to a constant value which means it needs 
longer period for the boundary system to form the Cooper pairs to complete the condensate. It again 
presents us the consistent picture on the influence of the nonlinearity in the electromagnetic fields on 
the scalar condensation. 
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